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^ , For a clean high-Tc superconductor, we analyze the Lawrence-Doniach 

^^1 free energy in a tilted magnetic field within the lowest Landau level (LLL) 

c/5 '_ approximation. The free energy maps onto that of a strictly c-axis field, but 

with a reduced interlayer coupling. We use this result to calculate the tilt 
modulus C44 of a vortex lattice and vortex liquid. The vortex contribution 
to C44 can be expressed in terms of the squared c-axis Josephson plasmon 
^ I frequency ojpi . The transverse component of the field has very little effect on 

O ' the position of the melting curve. 
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I. INTRODUCTION 



This paper is concerned with the tilt modulus C44 of the vortex system in a high-Tc or 
otherwise layered superconductor. C44 is an elastic constant which measures the free energy 
cost of applying a small transverse field in addition to a field applied parallel to the c axis. 
It is a relevant parameter in many physical processes, such as collective pinning |l|, the 
"peak effect," |l|J^, and the transition to the Bose glass state [^. Previous investigations 
have shown that C44 is very strongly dependent on the wave vector k of the transverse field 
0J^. In addition, it is finite in the flux liquid as well as the flux solid phase and is strongly 
affected by different kinds of disorder |I],^J31 

A number of authors have analyzed C44 in various approximations. Several groups 0,^,|^ 
have calculated C44(k) in a flux lattice as a function of the wave vector k of the transverse 
magnetic fleld. Other authors have considered C44 in the flux liquid state |0J^, and in the 
presence of disorder in both the solid and liquid phases |^,0 . 

In this paper, we calculate C44 for a layered superconductor at high fields. We describe 
the superconductor using a Lawrence-Doniach |10| free energy functional, and we evaluate 
fluctuations in the lowest Landau level (LLL) approximation, appropriate for strong c-axis 
magnetic flelds. This LLL approach accounts well for the position of the flux lattice melting 
curve in the magnetic-fleld/temperature plane [rT,r2|, as well as the value of the magnetiza- 
tion in both solid and liquid phases |TT| , p!3[| , both in Bi2Sr2CaCu208+5and the less anisotropic 
YBa2Cu307_a;. To obtain C44, we generalize the LLL approximation to apply to flelds with 
a flnite ab component. For a defect-free system in the LLL approximation, the tilted fleld 
free energy can be exactly mapped onto the usual free energy for a fleld B||c, but with a 
weaker interlayer coupling which depends on tilt angle. 

As a byproduct of this transformation, we can also analyze the behavior of the flux lattice 
melting curve in a magnetic fleld tilted at an angle to the c axis. We flnd that this curve is 
little affected by the presence of a transverse magnetic fleld component at flelds where the 
LLL approximation is applicable, consistent with available (but limited) experiment. 



II. LOWEST LANDAU LEVEL APPROXIMATION WITH TILTED FIELD 



At flxed external magnetic fleld H = HzZ + Hx^ the Lawrence-Doniach Gibbs free energy 
in a form which incorporates the fleld energy is |1^ 
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Here ipn is the order parameter in the n^^ layer, d is the distance between layers, q = — 2|e|, 
and t = h^ / {2mcd^) is the interlayer coupling energy. In a cuprate superconductor such 
as Bi2Sr2CaCu2084.5, d represents the repeat distance in the c direction. V_l is the xy 



component of the gradient operator, V is the system volume, and the magnetic induction 
B = V X A, where the vector potential A = (A_|_, A^). 

We choose the gauge A = —B^yk + B^yi and expand the ipn^ as a linear combination 

of ILL states, to obtain 7/;„(r) = { ^"'^^P \ Efc Cfc,„e^'=^-(^-'=^')'/2£2^ ^l^g^g ^ _ 27rp/L^ 



is the quantized momentum {p a positive integer), (i = J^qJI-rBz is the magnetic length, 
aniT) = [a(T) + 2t][l — Bz/Hc2(T)], and sample volume is V = L^LyLz- In terms of 
measurable quantities, a(T) + 2t = —hqHc2(T) / {2mabc) , where mab is the transverse effective 
mass and Hc2(T) is the upper critical field as a function of temperature T. Substituting this 
expansion into equation ([l|) and integrating over x and y yields 
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where N^ = (-Z^x/^)(v^/47r)^/^. This expresses the Gibbs free energy as a function of the 
expansion coefficients Ck^n for a tilted magnetic field. This expression differs from that for 
B||z []TT| in only two respects: (i) there is an extra phase factor in the interlayer coupling; 
and (ii) the strength of the interlayer coupling is renormalized by an exponentially decaying 
factor. 

The phase factor in eq. (^ can be eliminated by introducing a new set of coefficients 
Ck,n = Cfc,„e-*(2^^-'^^''=/*o)", in eq. (|). The product term becomes 4„Cfc,„_ie-*2^^-°'^''=/*o = 
Ckn^k,n-i- This transformation does not affect the free energy term which contains products 
of four coefficients. To see this, note that when a typical such product is transformed, 
it picks up a phase factor e*(2^^-^^'"/*o)('=-(^+5i)-('=+'?2)+{fc+gi+g2))^ ^hgre each {k + Qi) term 

is an integer modulo A^<j,, N^ being the number of flux lines in the sample. The terms 
in the exponential can be summed to yield zero, leaving this term unaltered. Similarly 

I |2 I ~ |2 

In terms of the new variables, the free energy is thus 
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where 



t' = te-(-^^'^^/*o)' (4) 



is the renormalized interlayer coupling. Thus, the tilted B-fleld free energy is identical in 
form to the free energy for a fleld B||z, but with a weaker interlayer coupling. The ground 
state solution for the redeflned amplitudes Ck,n are identical to those for the Ck,nS in a 
strictly longitudinal fleld but weaker B||z. However, the order parameter ipn picks up a B^ 
dependency. 
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corresponding, at low temperatures, to a tilted Abrikosov lattice. 

III. TILT MODULUS 

C44 is defined by 

where 6 is the angle between H and the c axis, and where Q, the Gibbs free energy appro- 
priate to an experiment at constant H and T, is given by 



g = -ksTln fl[ddl^dc,,ne-^/'-'' + V^-i^ ^ g^ + V^-^—^. (7) 
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Then, writing d/d9 -^ Hzd/dH^ for small H^, we find that the second term in eq. (|^ 
contributes a term C44 = Hi/ {An) ~ i?^/(47r) in the LLL regime, where the magnetization 
is small. This is the compressive contribution to the tilt modulus. 

To calculate the remaining (vortex-related) contribution to C44, which we denote C44, 
we first write G/ksT = Ti/T , where T = ^i ^^r^-Zp^ and 

T-i = N^}_^i -\Ck,n\ , T^, (Cfc^„Cfc,n_i + C.C.) 
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Then writing d"^ /d9^ = B'^d'^/dB'^, taking B^ ~ H^, and evaluating the derivatives from 
eqs. (^, with the result 

1 /9 
^^'' = [t) (^^^7^1^) ^^ - B^/H^^^^^(T.KnCk,n-^ + C.C.))t, (9) 

where we have dropped the tildes on the Cfc,n's because this expression is evaluated at B^ = 0. 
Eq. @ can be expressed in terms of measurable quantities using the identities a(T) + 2t = 
hqHc2(T)/2mabC and b = 2n h? {qh/niahc) , where k = Xab{Bz,T)/^abiBz,T) is the ratio of 
the ab penetration depth and coherence length. We also write t = ^^/(2mab(i^7^), where 7^ 
is the anisotropy parameter. The result is 

' * ^ k,n 

Finally, the total tilt modulus is 



^44 ~ ^ + CL. (11) 

In a triangular Abrikosov lattice, for example, only A^^^ = N^/N^ of the N^ coefficients 
in each layer are nonzero. The nonzero coefficients are C2pN.^,n = {'^1 PaY^"^-, C{2p+i)N^,n = 
i{2/ I3aY^'^i where p is an integer ranging from to Ny and /3^ = 1.159595... is the Abrikosov 
ratio |T^. Substituting these values into eq. (p!OD, using Ny/Ly = 2Nx:/{Lx^/S), and adding 
the compressive term, we obtain 

Cu = ^+ . ^\ , (gc2(T) - B,)B,. (12) 

Eq. (|12D can be compared to a previous estimate 0,1151 C^^'^ = (-B^/(47r)[l + 
l/(47rA^^(i?2,T)7^ns)], where Ug is the superfluid number density for a two-dimensional 
bose system related to the three-dimensional flux line system by a path integral formalism 
i. Writing Xab{B,,T) = kU{B,,T), 1/^^,(5,, T) ~ 27t{H,2{T) -B,)/^o, and approximat- 
ing Hs by riQ = Bz/^o, the number of bosons (i. e., flux lines) per unit area @J3]' ^^ obtain 
Cf/*^ ~ i?f/(47r) + [Hc2{T) - Bz)Bz]/{87TK'^'y'^) This resuh is in agreement with our own 
result, eq. (0), except for the factor (I/Pa) ~ 0.86, which arises from the nonuniformity of 
Hg in the Abrikosov lattice phase. 

Our C44 is closely connected to the so-called Josephson plasmon frequency ujpi, as cal- 
culated in the same LLL approximation [|16|. Ljpi{Bz,T) is the squared plasma frequency of 
the Josephson junction formed between adjacent ab layers, and is given by [Q 

u;liB,T) = iuj'J^r^[PANJi2N^N,)]{J2icl^c,,^., + cc))t, (13) 

k,n 

where u^^ = J[{Hc2{T) — B)cq]/[eoK'^'y^h(3A] is the mean-field Josephson plasmon fre- 
quency and eo is an interlayer dielectric constant. Combining eqs. (|10D and ([l3| ) gives 

CUB,,T) B^h' ^^^^ 
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Fig. 1 shows Cl^{Bz,T) for B^ = 2 tesla in clean Bi2Sr2CaCu2 08+<5, using this relation 
Upi as calculated in Ref. |jl6|. Like Upi, C44 has a discontinuity at the flux lattice melting 
transition, and remains flnite in the flux liquid state. Experiments on uj^i are consistent with 
this result, indicating that ujpi, and hence C44, have discontinuities at flux lattice melting in 
clean Bi2Sr2CaCu208+5 0>&- 



IV. FLUX LATTICE MELTING IN A TILTED MAGNETIC FIELD 

Eqs. (1^) and (^ have some striking implications for the LLL phase diagram of a clean 
layered superconductor. When B||c, this phase diagram depends on only two parameters. 



namely g = aH\lT^P'd/{bkBT) and r/ = t/la^l 131- ^^^ results show that, even with a 
transverse fleld B^-, the phase diagram still depends on only two parameters, except that r] 



is now replaced by r]' = r7exp[— (tt-B^c? /(2i?2$o)]. As found previously UTTHlq], this phase 



diagram contains a single first-order melting line separating a triangular vortex lattice from a 
vortex liquid. In the regime where the LLL approximation is adequate, our results show that 
this first-order line should persist in a tilted magnetic field (cf. Fig. 2). Furthermore, in most 
high-Tc materials, the line is shifted very little by a nonzero B^. (For Bi2Sr2CaCu208+5, 
t' ~ texp[-7.5 X lO-'^B^/B^].) 

We have found no experimental melting data in tilted magnetic fields at fields where the 
LLL approximation is applicable. In BiSr2Ca2Cu208+x, the low-field melting temperature 
has been reported |]19[ to depend only on the c component of H, which, though obtained 



in a very different regime, would be consistent with the result of our calculations. It would 
be of great interest to have further tests of the LLL predictions in the relevant high-field 
regime. 

V. SUMMARY 

In this paper, we have extended the LLL approximation for high-Tc superconductors 
to treat fields tilted at an angle to the layer perpendicular. The resulting free energy has 
exactly the same form as the usual case, except that the effective interlayer coupling is re- 
duced. For high-Tc materials, this reduction is small; hence, we predict that the flux lattice 
melting temperature will be little affected by the application of an oblique magnetic field in 
the range where the LLL approximation is valid, except possibly at angles nearly parallel 
to the layers. This prediction appears consistent with some existing (but low-field) experi- 
ments |jl9|. We also obtain an expression for the zero-wave- vector tilt modulus C44, in good 
agreement with previous estimates by other means P,P3|. Finally, the vortex contribution 
C44 is proportional to the squared Josephson plasmon frequency uj^i, as calculated in the 
same LLL approximation, and remains finite in the vortex liquid as well as the vortex solid 
phase. 
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Vortex contribution C44 



FIG. 1. 
Bi2Sr2CaCu208+5at B 



to the tilt modulus, plotted versus temperature for 



2T, as evaluated using eq. ([1 



and ujpi from Ref. 
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is discontinuous at the flux lattice melting temperature Tm{Bz) (indicated by arrow). Inset: 
enlargement of C44 near and above Tm- 

FIG. 2. Phase diagram of a clean high-Tc material in the LLL approximation, including a 
transverse magnetic field component B^. The parameters 1]' and g are defined in the text. 



Points and spline fit are given by Hu and MacDonald [12 
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